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ABSTRACT 

We present new warm Spitzer occultation photometry of WASP-26 at 3.6|J.m and 
4.5u\m along with new transit photometry taken in the g,r and i bands. We report 
the first detection of the occultation of WASP-26b, with occultation depths at 3.6/im 
and 4.5/oti of 0.00126 ± 0.00013 and 0.00149 ± 0.00016 corresponding to brightness 
temperatures of 1825 ±80K and 1725 ±89K, respectively. We find that the eccentricity 
of the orbit is consistent with a circular orbit at the la level (e = 0. 0028^° nn?9> 3cr 



upper limit e < 0.04). According to the activity- inversion relation of Knutson et al 
( |2010[ ), WASP-26b is predicted to host a thermal inversion. The brightness tempera- 
tures deduced from the eclipse depths are consistent with an isothermal atmosphere, 
although it is within the uncertainties that the planet may host a weak thermal inver- 
sion. The data are equally well fit by atmospheric models with or without a thermal 
inversion. We find that variation in activity of solar-like stars does not change enough 
over the tim e-scales of months or years to change the interpretation of the |Knutson| 



et al. (2010) activity-inversion relation, provided that the measured activity level is 
averaged over several nights. Further data are required to fully constrain the thermal 
structure of the atmosphere because the planet lies very close to the boundary between 
atmospheres with and without a thermal inversion. 

Key words: methods: data analysis - techniques: photometric - stars: individual: 
WASP-26. 



1 INTRODUCTION 

The first detection of thermal emission from an exoplanet 
was reported by Deming et al. ( 2005 1 and 



Charbonneau 



et al. (20051. The teams observed the secondary eclipse of 
HD209458 and TrES-1 using the Spitzer Space Telescope. 
Secondary eclipses of many other exoplanets have now been 
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observed (e.g. 


Machalek et al.||2008 Anderson et al.||2011b 


Todorov et al. 


2012 1. Through the spectrophotometry of this 



event, observed using Spitzer and ground based telescopes, 
we can build up the spectral energy distribution (SED) of 
the irradiated hemisphere (day side) of the planet. From the 
SED we can investigate the atmospheric properties of the 
day side of the planet. Secondary eclipse observations made 
with Spitzer have shown that some of these exoplanets have 



temperature inversions (Fortney et al.|2008 Knutson et al 
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2009||Madhusudhan fc Seager|2010) . Thermal inversions are 
thought to form when gases exist in the upper atmosphere 
of these exoplanets that are efficient absorbers of optical 
and ultraviolet light ( Fortney et al.||2008| ). This absorption 
of radiation causes the temperature of this region of the at- 
mosphere to increase. Gases that have been hypothesised 
to cause thermal inversions to form are titanium oxide and 



vanadium oxide ( Spiegel et al. 2009 ) and sulphur compounds 
( jZahnle et al.|2009[ ) 



WASP-26 b, discovered by |Smalley et al.| (|2010j with 
Super WASP ( |Pollacco et~aT1|2006| ), is a 1 Jupiter mass 
(lMj U p) planet in a 2.8 day orbit around a GO type star. 
WASP-26 also has a common proper motion companion 15" 



away (Smalley et al. 20101. Anderson et al. (2011al con- 



ducted an investigation using the Rossiter-McLaughlin ef- 
fect to determine the sky-projected spin-orbit angle of the 
system. However, their results were inconclusive. |Albrecht| 
et al. (20121 constrained the spin-orbit angle of the system 
to A = — 34l26°- in this paper we present new warm Spitzer 
and ground based photometry of WASP-26. 
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Figure 1. BUSCA and FTS data with best fitting models (from 
the top to bottom), g band data, r band data, i band data and 
the FTS (z filter) data. 



2 OBSERVATIONS 

We present Spitze r |Werner et al.||2004[ ) InfraRed Array 
Camera (IRAC) ([Fazio et al.||2004[) channel 1 (3.6um) 



and channel 2 (4.5(J.m) secondary eclipse (occultation) data 
taken on 2010 August 3 and 2010 September 7-8, respec- 
tively (PI: J H, Program ID 60003). The Spitzer data were 
acquired in full array mode (256 x 256). Also presented are 
new full transit data taken in the g, r and i bands (taken 
simultaneously) using the 2.2m telescope at the Calar Alto 
Astronomical Observatory with the Bonn University Simul- 
taneous CAmera (BUSCA) on 2010 August 20. BUSCA is 
a 4 channel CCD photometer with 4096 x 4096 pixels per 
CCD with a plate scale of 0.17 arc seconds per pixel. The 
BUSCA transit data were obtained using defocused photom- 



etry (Southworth et al.||2009c 20121 and BUSCA was used 



with a 256 x 1400 pixel window and 2x2 binning to reduce 
the read out time. Table [l] is a summary of the data that we 
have used in our analysis. 



3 DATA REDUCTION 



3.1 Transit Data Reduction 



We used an IDL implementation of DAOPHOT (Stetson 
19871 to perform synthetic aperture photometry on our 



BUSCA images, as in Southworth et al. ( 2009c I. Light curves 



were obtained in the g, r and i bands. In all three bands one 
comparison star was used. We used a target aperture radius 
of 24 pixels, a sky annulus of inner radius 70 pixels and an 
outer radius 100 pixels for the g, r and i bands. The wings 
of the PSF of the companion star do contaminate the target 
aperture but the contribution to the observed flux is negli- 
gible. Iterative outlier removal was used on the image values 
in the sky annulus to remove the effect of the light from the 
wings of the companion's PSF in the sky annulus. The light 
curves are shown in Figure [I] 



3.2 Spitzer Data Reduction 

The data reduction was conducted using the Image Reduc- 
tion and Analysis Facility (IRAF J] using the same method 
as | Anderson et aT] ( |2011b| ), described briefly below. 

We convert from MJy/sr to electrons using equation 
|TJ , where the gain, exposure time and flux conversion factor 
were taken from the image headers. 



Factor = 



Gain x Exposure Time 
Flux Conversion Factor 



(1) 



Aperture photometry was then conducted using the 
PHOT procedure in IRAF, using 21 aperture radii in the 
range 1.5-6 pixels and with a sky annulus of inner radius 8 
pixels and outer radius 16 pixels. It was found that the stel- 
lar companion to WASP-26 and a bad column in channel 
2 data were both inside the sky annulus. However, an iter- 
ative 3-sigma clipping was conducted which excludes those 
pixels. The error on the photometry was calculated from 
the photon statistics and the read out noise of the IRAC de- 
tectors. The readout noise values were taken from the IDL 
program SNIRAC_warm.proQ the values for channel 1 and 
2 are 9.87 and 9.4 electrons, respectively. The position of the 
target was measured by fitting a 1-dimensional Gaussian to 
the marginal distributions of flux on x and y image axes. 
For each data set the times of mid-exposure were converted 



to BJDtdb (Eastman et al. 20101 and for the occultation 



data the light travel time across the system (~40s) was ac- 
counted for. The light travel time across the system was 
calculated using the semi-major axis from the output of our 
initial MCMC (see below for details of this run) and this 
time was subtracted from all the Spitzer times. 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

2 ssc.spitzer.caltech.edu/warmmission/propkit 
/ som I snirac- warm.pro 
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Table 1. Summary of data used in this analysis. 



Observation 



Dates 



Publication 



SuperWASP Lightcurves (400 - 700 nm filter) 

16 RV spectra from CORALIE 
(1.2m Swiss Telescope, La Silla, Chile) 
30 RV spectra from HARPS 

(HARPS Spectrograph, ESO 3.6m telescope, La Silla, Chile) 
Full transit (Pan-STARRS-z filter) 

(2.0m Faulkcs Telescope South, Siding Spring, Australia) 

Occultation (3.6/im) 

{Spitzer channel 1) 

Occultation (4.5/^m) 

{Spitzer channel 2) 

Full transit {g, r and i band) 

(Calar Alto Astronomical Observatory with BUSCA, Almeria, Spain) 



2008 June 30 - 2008 November 17 

2009 June 28 - 2009 November 17 

2009 June 19 - 2009 August 22 

2010 September 12 

2009 November 18 

2010 August 3 

2010 August 7 - 2010 August 8 
2010 August 20 



Smalley et al. 1 2010 I 



Smalley et al. 



(20101 



Anderson et al 



(2011a) 



Smalley et al. 



(2010) 



This Paper 
This Paper 
This paper 



4 ANALYSIS 
4.1 MCMC 

We explored the parameter space using a Markov chain 
Monte Carlo (MCMC) algorithm (|Collier Cameron et al.| 



Table 2. Proposal parameters of the model used in our MCMC 
analysis 



2007||Pollacco et al.|2008||Enoch et al.|2010[ ). The input pa- 
rameters for the star that were used in the MCMC analysis 
areTeff = 59 50 ±100 and [Fe/H] = -0.02±0.09 ( |Anderson| 



et al. 2011al. Stellar density, which is directly constrained 



by the transit light curve and the spectroscopic orbit ( Sea- 
ger fc Mallen-Ornelas|[2"003 1 and the eccentricity of the or- 



bit, is calculated from the proposal parameter values. This 
is input, together with the latest values of T e g and [Fe/H] 
(which are controlled by Gaussian priors) into the empirical 



mass calibration of Enoch et al. ( 2010 1 to obtain an estimate 



of the stellar mass, M*. At each step in the MCMC proce- 
dure, each proposal parameter is perturbed from its previ- 
ous value by a small, random amount. From the proposal 
parameters, model light and RV curves are generated and 
X 2 is calculated from their comparison with the data. A step 
is accepted if x 2 (our merit function) is lower than for the 
previous step, and a step with higher \ 2 i s accepted with 
probability exp(— Ax 2 )- In this way, the parameter space 
around the optimum solution is thoroughly explored. The 
value and uncertainty for each parameter are taken as the 
median and central 68.3 per cent confidence interval of the 
parameter's marginalised posterior probability distribution, 
respectively ( Ford|2006 ). The median closely approximates 
the x 2 minimum for symmetric posteriors such as ours, and 
is more robust to noise in the case of flat minima. Table 
[2] show the proposal parameters of the MCMC. We did an 
initial run which included all the transit photometry, in- 
cluding WASP photometry, to get a good estimate of the 
epoch of mid-transit. This value along with its uncertainty 
were used as a Bayesian prior in subsequent MCMC runs 
which used all the photometry, including Spitzer, but ex- 
cluding the WASP photometry (to reduce computing time). 
The transit model used in the analysis was the small planet 



approximation of Mandel & Agol ( 2002 1 with 4-parameter 
limb darkening coefficients taken from Claret (2004). The 
limb darkening coefficients were determined using an initial 
interpolation in logg* and [Fe/H] and an interpolation in 



p 

AF 
b 

Ki 
T e ff 

-y/ecosoj 
y/esmu) 
AF 3 . 6 
AF4.5 



Time of mid transit 

Period of planet 

Depth of transit 

Transit duration 

Impact parameter 

Stellar radial reflex velocity 

Effective temperature of the star 

Metallicity of the star 

e=eccentricity, u = argument of periastron 

Depth of secondary eclipse at 3.6^m 
Depth of secondary eclipse at 4.5^m 



T c ft at each MCMC step. The limb darkening parameters 
used for the best-fit lightcurves are given in Table [3] For the 
secondary eclipse we approximated the star and planet as 
two uniform discs of constant surface brightness. We fixed 
the projected spin-orbit angle to the value A = in our fit 
since the HARPS data covering the transit are negligibly 
affected by the R-M effect. The fit to the optical lightcurves 
(Figure [TJ shows that there is some red noise present in 
the g and i band lightcurves. We have accounted for the 
small additional uncertainty due to this noise in our quoted 
parameter standard errors rather than trying to find an ar- 
bitrary model that would improve the fit. 

We checked for any correlations in our proposal param- 
eters and only found the correlation between transit depth, 
width and impact parameter often seen in ground-based 
lightcurves. These correlations are caused by the blurring of 
the second and third contact points due to limb darkening in 
the optical lightcurves. These correlations do not affect our 
secondary eclipse depth measurements. These correlations 
are shown in Figure [2] We also checked that our chain had 
converged, both by visual inspection and using the Gelmen- 
Rubin (G-R) statistic ( |Gelman et al.|2003| |Ford|2006[ ). 
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0.74 0.76 0.78 0.80 0.82 0.84 0.86 
Impact Parameter 




0.74 0.76 0.78 0.80 0.82 0.84 
Impact Parameter 



0.034 0.035 0.036 0.037 
Transit Duration 





0.74 0.76 0.78 0.80 0.82 0.84 0.E 
Impact Parameter 



0.74 0.76 0.78 0.80 0.82 0.84 0.86 
Impact Parameter 



Figure 2. Correlation plots for selected proposal parameters from our MCMC analysis. For clarity we have only plotted a random 2% 
of the chain values. 



Table 3. Limb darkening coefficients 



Light Curve 


"i 


a 2 


03 


a 4 


FTS 


0.655 


-0.352 


0.645 


-0.329 


BUSCA (g band) 


0.433 


0.208 


0.496 


-0.300 


BUSCA (r band) 


0.555 


0.028 


0.445 


-0.278 


BUSCA (i band) 


0.641 


-0.267 


0.640 


-0.338 



4.2 Trend Functions and Aperture size 

Figure [3] shows an example of the 3. 6 urn light curve pro- 
duced by the photometry in IRAF. There is a steep increase 
in the measured flux during the first part of the observation. 
This occurs because the telescope has slewed from its old po- 
sition to its new position and is adjusting to a new equilib- 
rium. We exclude the data that precedes HJD=2455447.37, 
to remove the major part of the initial ramp. It can be seen 
that there is a clear periodic trend in the data. This is due 
to the variation in the position of the target on the detec- 
tor caused by flexure of the instrument as an electric heater 
is turned off and on0 The IRAC detectors are known to 



exhibit inhomogeneous intrapixel sensitivity (e.g. Knutson 
et al.||2008[ ), which means that different parts of the detec- 
tor are more or less sensitive than others. This, along with 
the PSF movement, results in the measured flux varying de- 
pending on the position of the PSF on the detector. Also, 
when small apertures are used pixelation occurs due to the 



ssc.spitzer. caltech. edu/warmmission /news/21oct201 Omemo . pdf 



£ 1.22x10 5 k + . 




5446.90 5446.95 5447.00 5447.05 5447.10 5447.15 
HJD-2450000.5 



Figure 3. The raw light curve of the 3.6u.m Spitzer data ex- 
tracted using and aperture of 2.4 pixels. 



under-sampling of the PSF of the target (Anderson et al 



2011b I. These systematics will be accounted for in the trend 
functions as described below. Figure [4] shows an example of 
the 4.5um data which is less affected by these systematics 
even though (as it can be seen from Figure [5]l the radial 
motion of the PSF is greater at 4.5u.m than at 3.6um. 

The general form of the trend functions that were used 
in our analysis is 



Af = a + a x Ax + a y Ay + a xy AxAy + a xx Ax 2 
+a yy Ay + at At, 



(2) 
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Figure 5. The top, middle and bottom plots of the each column show the distance of the PSF from the nearest pixel centre in x, y and 
radially in each of the measured wavelengths respectively. 
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+ + 




5411.05 5411.10 5411.15 5411,20 5411.25 5411.30 
HJD-2450000.5 

Figure 4. The raw light curve of the 4.5|J.m Spitzer data ex- 
tracted using an aperture of 2.4 pixels. 



parameters we calculate the residuals from the model and 
then calculate the coefficients of the detrending model using 
singular value decomposition applied to the entire data set. 
Initially, a linear-in-time and quadratic-in-space trend func- 
tion was used on all 21 apertures to fit the secondary eclipse 
data. The RMS of the residuals was used to determine the 
optimal aperture size. Once this was determined, combina- 
tions of no, linear and quadratic trend functions in time and 
space were used on the best aperture to determine the best 
fitting trend function. 

Initially this decorrelation was conducted using the po- 
sitions measured by the 1-dimensional Gaussian fit to the 
target. We also attempted to remove the trends in the data 
by decorrelating against the radial position (radial distance 
from the centre of the nearest pixel) instead of the x and y 
positions independently. The general trend function for the 
radial decorrelation is, 

A/ = b + bir + b 2 r 2 + b t At, (3) 



where A/ = / — / is the stellar flux relative to its weighted 
mean, Ax — x — x and Ay = y — y are the coordinates of the 
point spread function of the target centre relative to their 
weighted means, At is the time since the beginning of the 
observation, and an, a x , a y , a xx , a yy and a t are coefficients 
which are free parameters in the MCMC analysis (Ander- 



son et al. 2011b I. For each set of trial lightcurve model 



where bo, 6i, 62, bt are free parameters in the MCMC analysis 
and r is the radial distance from the centre of the nearest 
pixel centre. A third method that was attempted was to use 
target positions in the trend functions measured by fitting 
a two dimensional circular Gaussian of fixed full width half 
maximum (1.39 pixels in channel 1 and 1.41 pixels in channel 
2) to a small region of the images containing the target. 
To determine which trend function gave better results 
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we used the Bayesian Information Criterion (BIC) (Schwarz 
19781), 



BIC = x + fclnO) 



(4) 



where k is the number of free parameters and n is the num- 
ber of data points. This method of determining how com- 
plicated a model to use only accepts a higher order trend 
function if the fit improves \ 2 by ln(n) or better for each 
additional free parameter. 

Using the RMS of the residuals it was found that the 
best aperture to use was 2.4 pixels in both channels. The 
system parameters are negligibly affected by the choice of 
aperture radius around this value. It was also found that the 
RMS of the residuals to the channel 1 data was marginally 
lower when using the position measurements measured by 
the 2D circular Gaussian method as opposed to ID Gaussian 
position measurements (0.002995 compared to 0.003054). 
The channel 2 data gave consistent RMS no matter the po- 
sition measurement used. The system parameters were con- 
sistent no matter which position measurement system were 
used. The results shown in Figure [6] and Table [| are those 
using the 2D circular Gaussian method, extracted from the 
2.4 pixel aperture and trend functions as described below. 
We found that the radial decorrelation gave a worse fit to our 
data compared to that of x and y decorrelation (x 2 worse 
by ~ 3000 at 3.6um and ~ 400 at 4.5um). 

Using equation Q it was found that the quadratic-in- 
space with no time trend function gave the best fit to the 
data in channel 1 and that the linear-in-space with no time 
trend function gave the best fit to the data in channel 2. It 
was found that the addition of the quadratic term for the 
spacial decorrelation improved our BIC by ~ 200 in channel 
1 and less than ~ 10 for more complicated models in both 
channels. We also detrended our data based only on the out- 
of-eclipse points to see if this affected our measured eclipse 
depths. It was found that the eclipse depths were consistent 
with our previous decorrelation. 



5 RESULTS AND DISCUSSION 

5.1 Eclipse Depths and Brightness Temperatures 

We find that the eclipse depths at 3.6u.m and 4.5(im are 
0.00126±0.00013 and 0.00149±0. 00016, respectively. These 
eclipse depths correspond to brightness temperatures of 
1825±80K and 1725±89K. To find these blackbody temper- 
atures the expected flux ratios were calculated using Planck 
functions at different temperatures for the planet and syn- 



thetic spectra from stellar models (Kurucz 1991 1 for the star 



These flux ratios were then integrated over the Spitzer band 
passes to calculate the expected measured flux ratio. The 
temperatures above correspond to the best fitting Planck 
function temperature to the individual eclipse depths. The 
errors were calculated using a simple Monte Carlo method. 
These temperatures suggest that, on average, the emission at 
mid-infrared wavelengths from the irradiated hemisphere of 
WASP-26b is consistent with the spectrum of an isothermal 
atmosphere, with the possibility of a weak thermal inversion 
within the uncertainties on the brightness temperatures. 



5.2 Atmospheric Analysis 

We model the day-side emergent spectrum of the hot Jupiter 
WASP-26b using the atmospheric modeling and retrieval 
technique of |Madhusudhan fc Seager] ( |200l)} [2010] ). The 
model computes line-by-line radiative transfer in a plane- 
parallel atmosphere in local thermodynamic equilibrium, 
and assumes hydrostatic equilibrium and global energy bal- 
ance. The pressure-temperature (P-T) profile of the atmo- 
sphere and the chemical composition, i.e. the sources of 
molecular line opacity, are input parameters to the model. 
The model atmosphere includes the major sources of opacity 
expected in hot, hydrogen-dominated atmospheres, namely, 
molecular absorption due to H2O, CO, CH4, and CO2, and 
continuum opacity due to H2-H2 collision-induced absorp- 



tion (CIA). Our molecular line-lists are discussed in Mad- 
husudhan fc Seager| | |2009| | and |Smith et al.| ( |2012| |. Given a 
photometric or spectral dataset of thermal emission from the 
planet, we explore the space of atmospheric chemical compo- 
sition and temperature structure to determine the regions in 
model space that explain, or are excluded by, the data (e.g. 
Madhusudhan et al.|2011 1. In the present case, however, the 
number of available data points (N = 2) are far below the 
number of model parameters (N — 10), implying that a 
unique model fit to the data is not feasible. Consequently, 
we nominally fixed the chemical composition of the models 
to that obtained with solar elemental abundances in ther- 
mochemical equilibrium (e.g. Burrows & Sharp 1999; Mad- 



husudhan 2012) for a given thermal profile, and explored the 



space of thermal profiles, with and without thermal inver- 
sions, that might explain the data. 

Figure [7] shows the 3.6um and 4.5u.m data along with 
model spectra of atmospheres with and without a thermal 
inversion, and a blackbody model. All three models shown 
allow for very efficient day-night redistribution. We find that 
both our planet-star flux ratios can be explained by a plan- 
etary blackbody at around 1750 K. Consequently, the data 
are consistent with an isothermal atmosphere. However, an 
isothermal temperature profile may be unphysical in radia- 



tively efficient atmospheres at low optical depth (e.g. Hansen 
2008). A temperature profile with a non-zero thermal gra- 
dient, with or without a thermal inversion, may be more 
plausible. As shown in Fig. [7j the two data are fit almost 
equally well by models with and without a thermal inversion, 
as shown by the red and green models, respectively. Further 
occultation depths measured at different wavelengths are re- 
quired to break the degeneracies between the models and to 
determine the true nature of the atmosphere. It can be seen 
in Fig. [7]that there are some differences between the models 
with and without a thermal inversion at 1.25um (J band), 
1.65u.m (H band) and 2.2u.m (K band). These wavelengths 
are accessible from the ground, so with measurements of 
the occultation depth at these wavelengths it may be pos- 
sible to break the degeneracies between these models. Hub- 
ble Space Telescope WFC3 observations covering the wave- 
length range 1 — 1.7u.m can also be used to detect spectral 
features due to water either in emission or absorption, and 
so distinguish between models with and without a thermal 
inversion ( |Madhusudhan||2012| [Swain et al.||2012[ ). We em- 
phasize that we have only presented two possible models 
here that represent the average properties of the irradiated 
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Figure 6. (Left): The raw light curves with the trend functions, the upper points are the channel 2 data and the lower points are the 
channel 1 data, the solid lines are the trend functions for each data set. (Middle) Binned light curves with trend models. (Right) The 
binned light curve with trend function removed and best fitting eclipse models (solid lines). The secondary eclipse can clearly be seen in 
both channels. 
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Figure 7. Spectral energy distribution of WASP-26b relative to 
that of its host star. The blue circles with error bars are our best- 
fitting occultation depths. The green line is a model-atmosphere 
spectrum, based on a model which assumes solar abundances in 
thcrmochemical equilibrium and lacks a temperature inversion, 
and the dark red line is a model with a temperature inversion. The 
band-integrated model fluxes are indicated with circles of the cor- 
responding colours. The dashed black line shows a planetary black 
body model with a temperature of 1750K. Inset: temperature- 
pressure profiles for our models. 



hemisphere of WASP-26b. With additional data other pa- 
rameters of the models such as composition can be explored. 



Knutson et al. (2010 1 (hereafter K10) presented results 



which suggest that planets without thermal inversions orbit 
active stars, and those with inversions orbit inactive stars. 
This may be due to photodissociation of the opacity source 
in the upper atmosphere of the planet by the UV flux from 
the active stars (K10). It is known that solar-like stars have 
activity cycles on time scales of approximately 10 years. The 



Duncan et al. ( 1991 I catalogue of Shk activity measure- 



ments taken at the Mount Wilson Observatory was used to 
examine to what extent the activity of a star changes on 
short time scales (order of months) and long time scales (or- 
der of years). The aim was to determine if the variability in 
activity of the stars in the K10 sample was such that, in the 
time between the occultation observation and the measure- 
ment of \ogR' HK , the activity of the star can change enough 
to affect the interpretation of this activity-inversion relation. 



Recently Montalto et al. (20121 showed that the activity of 
WASP-3 changed from \ogR HK = —4.95 (less active) to 
\ogR' HK = -4.8 (more active) between 2007 and 2010. It 



has been shown by Menou & Rauscher ( 2009 ) that the time 



scale for models of hot Jupiter atmospheres to go from their 
initial conditions to a statistical steady state was ~ 20 days. 
This suggests the time scale of hot Jupiter atmosphere vari- 
ability is much shorter than the time scale of stellar activity 
variability. More detailed modelling and additional observa- 
tions are required to better understand whether variations 
in the UV irradiation can produce observable changes in the 
eclipse depths for planets near the boundary between atmo- 
spheres with and without strong thermal inversions. 
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Table 4. System parameters from our MCMC analysis 



Parameter 



Symbol (unit) Value 



Orbital period 

Epoch of mid-transit (BJD, TDB) 

Transit duration (from first to fourth contact) 

Duration of transit ingress duration of transit egress 

Planet-to-star area ratio 

Impact parameter 

Orbital inclination 

Semi- amplitude of the stellar reflex velocity 
Centre-of-mass velocity 

Argument of periastron 



Orbital eccentricity 
Phase of mid-occultation 
Occultation duration 
Duration of occultation ingress 

Star mass 

Star radius 

Star surface gravity 

Star density 

Star effective temperature 
Star metallicity 

Planet mass 

Planet radius 

Planet surface gravity 

Planet density 

Semi-major axis 

Occultation depth at 3.6/^m 

Occultation depth at 4.5/^m 

Planet equilibrium temperature (full redistribution)* 
Planet equilibrium temperature (day-side redistribution)* 
Planet equilibrium temperature (instant reradiation)* 



duration of occultation egress 



P(d) 
T c 

T U (d) 
T 12 « T 34 (d) 
AF = Rl/Kl 
b 

t(°) 

Ki (km s- 1 ) 
7 (km s _1 ) 

w (°) 
ecosa; 
es'mui 
e 

^mid — occultation 

T 58 (d) 

T 56 « T 78 (d) 

M, (M Q ) 
R* (Ha) 
logs* (cgs) 

p* (pa) 
T e g (K) 
[Fe/H] 

M P (Af Jup ) 
flp (Rjup) 
loggp (cgs) 

Pp (pj) 
a (AU) 
AF3.6 
AF4.5 

5p,A=o,f=i (K) 

Tp,A=0,f=2 (K) 

(K) 



2.756611 ±0.000008 

2455424.10899 ±0.00012 

0.097 ±0.002 

0.024 ±0.002 

0.0103 ±0.0003 

0.82 ±0.02 

82.9 ±0.4 

0.138±0.002 
8.4593 ±0.0001 

+200 



"90-20 
-0.0004 ±0.0007 



-0.0011 



P,A=0,f= 



+0.0023 
-0.0110 
n nmsQ +0.00965 
0.00283 _ 002 2i 

0.4998 ±0.0005 

0.097 ±0.002 

0.024±0.002 

1.10±0.03 
1.29±0.05 

4.26 ±0.03 
0.52±0.06 
6000±100 
-0.02±0.09 

1.03±0.02 

1.27 ±0.07 
3.16 ±0.04 
0.50 ±0.08 
0.0398 ±0.0003 
0.00126 ±0.00013 
0.00149 ± 0.00016 
1623 ±43 

1930 ± 51 
2074 ± 55 



* where A is the albedo, f=l is defined as full redistri bution, 

f=2 is day-side redistribution, and f=| is instant reradiation as in Smith et al. (2011 1 



We converted the Shk measurements in Duncan et al 



( 199l[) to \ogR' HK using the method described by Noyes 



et al.l ( 1984 1. A look-up table based on \ogR' HK and B-V 



colour for the stars in the Duncan et al. ( 1991 1 catalogue was 



then constructed. Using this table, the within-season varia- 
tion of log R' H K of the stars was used as a measure of the 
short term variability in log_R^ K and the season-to-season 
variation in log R'hk as a measure of the long term variation 
in \ogR' HK . This look up table was then used to estimate 
the variation in \ogR' HK for the stars of K10 based on their 
B-V colour. It was found that the short term variability 
was always ^ 0.02 dex and the long term variability was be- 
tween 0.02 and 0.06 dex. This suggests that the variation in 
\ogR' HK is not large enough on either short nor long term 
time scales to change the interpretation of K10. However, 
this may blur the boundary between the two classes of plan- 
ets. The error bar shown in Figure [8] is the typical change 



in activity, assuming the spectra are measured over several 
nights. It is possible for stars to vary by much more than 
this amount over their rotation period (e.g. |Dumusque et al.| 
|2012[ |. This short time scale variation will move the star on 
the diagram but this may not reflect changes in UV irradi- 
ation. The value of log R' HK = -4.98 for WASP-26 used in 



(2011a) 



this analysis is taken from Anderson et al. 

We compiled updated values of R p /R* and the sec- 
ondary eclipse depths for the stars in the K10 sample. Figure 
5 of K10 was then replotted, this is shown in Figure [8] We 
include on this plot WASP-26b. As can be seen from Fig- 
ure |8j it seems to lie very close to the boundary between 
the two classes. Using the convention as in [Anderson et al.| 
( |2011c[ ) the abscissa value for WASP-26b is < = -0.020 ± 
— 0.023%um _1 , where £ is the gradient of the measurements 
at 3.6um and 4.5(J.m, i.e. AFs.e — AF4.5/(— 0.9um), minus 
the gradient of the blackbody that is the best-fit to the 
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two measurements. The theory behind this is that at 4.5u.m 
there are opacity sources that are not present at 3.6(im (CO 
and H 2 0) ( |Madhusudhan fc Seager||2010| ). The 4.5um data 
probes a higher region of the atmosphere compared to the 
3.6um data. This suggest that if the brightness temperature 
at 4.5u.m is greater than that at 3.6u.m then there is likely 
to be a thermal inversion in the atmosphere. 



5.4 Ecentrictiy 

From secondary eclipse measurements it is also possible to 
constrain the eccentricity of the orbit from timing of the 
secondary eclipse relative to transit. We find that the ec- 
centricity of the orbit is small (e — 0.00281q'qo22)i which is 
consistent with a circular orbit at the la level. We find a 3cr 
upper limit on the eccentricity of the planet's orbit of 0.0399 



which is similar to Anderson et al. (2011a I 3er upper limit 
of 0.048. 



6 CONCLUSION 

In this paper we present new warm Spitzer photometry 
of WASP-26 at 3.6um and 4.5u.m along with new tran- 
sit photometry taken in the g,r and i bands. We report 
the first detection of the occultation of WASP-26b with 
eclipse depths at 3.6u.m and 4.5u.m of 0.00126 ±0.00013 and 
0.00149 ± 0.00016 respectively which correspond to bright- 
ness temperatures of 1825 ±80K and 1725 ±89K. Our anal- 
ysis shows that the atmosphere of WASP-26b is consistent 
with an isothermal atmosphere with the possibility of a weak 
thermal inversion (within the uncertainties on the bright- 
ness temperatures). If the K10 activity- inversion relation 
holds for WASP-26b, then we would expect it to host a 
thermal inversion. More secondary eclipse data at different 
wavelengths, particularly near-IR secondary eclipse depths 
near the peak of the planet's SED, will be able to better 
constrain the true nature of the atmosphere of WASP-26b. 
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Figure 8. Activity-inversion plot for the stars in |Knutson et ah] | |2010| l. Points on the left of the dotted line (triangles) are non-inverted 
planets around active stars and those on the right of the dotted line (squares) are inverted planets around inactive stars. The point on 



the left hand side of the plot shows the typical change in ' 
blue circle is WASP-26. The stars are listed in table [5] 



' y R'„ K (season-to-season), assuming measurements over several nights. The 



Table 5. Stars in Figure [8] 



Star 



Log(R' K ) 



Q value 



Rp/R* 



Eclipse depth in channel 1 Eclipse depth in channel 2 



Carter fc Winn jgOlO i 
Southworth (2011) 
South-worth ( 20(55 1 



HD189733 

TRES-3 

TRES-1 

WASP-4 

XO-2 

TRES-2 

XO-1 

HAT-P-1 

HD209458 

TRES-4 

COROT-1 

WASP-1 

WASP-2 

WASP-18 

HAT-P-7 

WASP-26 



-4.501* 
-4.549* 
-4.738* 
-4.865* 
-4.988* 
-4.949* 
-4.958* 
-4.984* 
-4.970* 
-5.104* 
-5.312* 
-5.114* 
-5.054* 
-5.430* 
-5.018* 
-4.98** 



-0.1707 
-0.0721 
-0.0682 
-0.0961 
-0.0271 
0.0447 
-0.0135 
0.0178 
0.0770 
-0.0301 
-0.0343 
0.0599 
0.0251 
0.0332 
0.0481 
-0.0200 



± 0.0271 
± 0.0715 
± 0.0204 
± 0.0457 
± 0.0292 
± 0.0354 
± 0.0127 
± 0.0260 
± 0.0194 
± 0.0216 
± 0.0660 
± 0.0259 
± 0.0432 
± 0.0176 
± 0.0309 
± 0.0229 
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